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Abstract 

BSAURUS is a software package for the inclusive reconstruction of B-hadrons in 
Z° bb events taken by the DELPHI detector. The BSAURUS goal is to recon- 
struct B-decays, by making use of as many properties of b-jets as possible, with 
high efficiency and good purity. This is achieved by exploiting the capabilities of 
the DELPHI detector to their extreme, applying wherever possible physics knowl- 
edge about B production and decays and combining different information sources 
with modern tools- mainly artificial neural networks. This note provides a reference 
of how BSAURUS outputs are formed, how to access them within the DELPHI 
framework, and the physics performance one can expect. 



1 Introduction 



BSAURUS is a package for inclusive B reconstruction of decays into b quark- ant iquark 
pairs taken by the DELPHI detector. Due to the large mass of B hadrons there are many 
thousand decay channels, all with small branching ratio. This is a severe limitation for 
studying B physics. BSAURUS tries to reconstruct inclusively as many properties of 
b-jets as possible with high efficiency and yet with good purity also. This is achieved 
by exploiting the capabilities of the DELPHI detector to their extreme, applying wher- 
ever possible physics knowledge about B production and decays and combining different 
information sources with modern tools- mainly artificial neural networks. 

1.1 Physics Motivation 

High energy b-jets contain B+,B'', Bg mesons or B-baryons in an approximate mixture of 
40:40:10:10. To measure hfetimes one needs to be able to reconstruct the decay length, 
i.e. the distance between primary and secondary vertex, as well as the B-hadron energy. 
The mean lifetime of B-hadrons have been measured with good accuracy and together 
with the scmilcptonic branching ratio it can be used to determine the Cabibbo-Kobayashi- 
Maskawa matrix element Vcb- The lifetimes of the individual B hadrons are predicted by 
the Heavy Quark Effective Theory (HQET) to be only marginally different from the mean 
lifetime. In order to test HQET, one needs to be able to distinguish between the four B 
hadron species. 

The fraction of the available energy that B-hadrons attain after the fragmentation 
process is not directly calculable in perturbative QCD, and has to be phenomenologically 
modelled with fragmentation functions. Measuring the b-quark fragmentation function, 
ideally for the different B-types separately, is fundamental to the accurate simulation of 
B-physics. 

The decay into b-quarks exhibits a forward-backward asymmetry between the b 
and b-jet. A precise measurement of A°p^ constitutes an important test of the Standard 
Model. For an experimental study one needs to tag the production fiavour (i.e. decide 
whether a given jet is from a b-quark or b antiquark). Neutral B mesons can oscillate into 
their antiparticlc (and back) before they decay. The time dependence of Bj oscillations has 
been observed, whereas the oscillations of the Bg meson are too fast to be experimentally 
resolved at present . For experimental studies of B°-oscillations one needs to know the 
production fiavour and the decay fiavour. 

The weakly decaying B hadrons are the lowest mass, i.e. the ground states, of hq 
systems with quantum numbers J^'^ = ^. For each of these systems a full spectroscopy 
of excited states exists, e.g. B* vector mesons and orbital excitations which are predicted 
to appear as two doublets for each unit of orbital angular momentum L, one being broad 
decaying in S-wave, and the other narrow since at least D-waves are necessary in the 
decay. Signals of L = 1 B-mesons are observed with a relatively large rate (about 20- 
30%), but the composition of the broad peak observed is still unclear. Radial excitations 
are also predicted in the quark model. In B-baryon spectroscopy, in addition to the Ab, 
the states S^, and Qi, are expected to decay via weak interactions with an appreciable 
lifetime, whereas the Eb and should decay into AbTr via strong interactions. For 
testing B-spectroscopy one needs the reconstruction of B-hadron four-momenta and the 
identification of protons, pions and kaons originating from the fragmentation process. 
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BSAURUS is designed primarily to allow analyses in all these areas. The algorithms 
used rely on a satisfactory performance from all of the detector components in the DELPHI 
barrel in conjunction with optimised reconstruction software. 

1.2 Detector-Specific Reconstruction: The Basis of BSAURUS 

Charged particle reconstruction, especially the correct assignment of silicon vertex detec- 
tor points to tracks reconstructed in the TPC, the Inner Detector and the Outer Detector, 
has been optimised in the years 1995/1996 by better alignment procedures and adding 
completely new algorithms including e.g. inside-out tracking, a complete event ambigu- 
ity resolution, and a Kalman filter track fit taking into account multiple scattering and 
energy loss. All LEP I data taken from 1992-1995 have been reprocessed with these new 
algorithms, and this has led to huge reconstruction improvements, especially in dense 
jets. As an example, the signal of exclusively reconstructed D* mesons into D^tt^ with 

— > KTTTTn was increased by a factor 2.4. In addition, optimised photon conversion 
algorithms [|l| , reconstruction of secondary hadronic interactions in detector material and 
of H~ decays, and a search for additional low-energy tracks in the vertex detector 
help to give a clean reconstruction of complicated events. It should be noted that since 
BSAURUS has been developed and tuned on these reprocessed data the result of applying 
it to any other data set cannot be guaranteed. 

Particle identification is an important ingredient for many of the algorithms applied. 
To achieve an optimal performance, electron identification [0 is performed using a neural 
network combining spatial and energy information from the HPC calorimeter, tracking in- 
formation, dE/dx measurements of the TPC, and searches for kinks in the track in known 
layers of material and for photons radiated tangentially in them. In DELPHI, hadron 
identification information is gathered from the Ring Imaging Cerenkov Counters(RICH) 
with both gas and liquid radiator, and from energy loss measurements in the TPC and the 
silicon vertex detector. To optimise hadron classification at various momenta, detailed 
information from all of these detectors has been combined using neural networks in the 
program package MACRIB ||^. 

1.3 Overview Of BSAURUS Outputs 

BSAURUS provides a host of information relating to the properties of B-hadron produc- 
tion, fragmentation and decay within an event hemisphere. This information is available 
via output variables stored in COMMON blocks, described in Appendix C. 
Among the most important outputs for B-physics analysis topics are: 

• estimates of the B-hadron energy(see Section |3.3| ) and decay length(see Sec- 
tion |3j), 

• the TrackNet (see Section ^ that discriminates between tracks originating from 
the primary and a secondary vertex, 

• the BDnet (see Section |^ that discriminates between tracks originating from the 
B and the cascade B^D vertex, 

• neural networks (see Section ^ that provide a probability that a particular B- 
hadron type (i.e. B"'",B°,Bs or B-baryon) was present in the hemisphere. 
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• neural networks (see Section ^ that tag the charge or flavour of the b-quark in 
the B-hadron type at both the production and decay time. 

2 BSAURUS Definitions 

2.1 Multihadronic Events 

There is no event selection applied inside BSAURUS i.e. it is left to the user to decide 
on which events to call BSAURUS. However, for all development and tuning work and 
for all results presented in this note, events in real data and Monte Carlo had to pass 
a multihadronic event selection. The criteria used are from the DELPHI pilot record 
hadron selection bits: 

• team 4 selection, > 5 charged particles (bit 1) 

• team 4 selection, charged energy > 0.12Elep (bit 2) 

• Eforward,Ebackward > 0.03ElEP (bit 12). 

2.2 Standard Particle Selection 

Charged particles used by the BSAURUS algorithms must pass the following track selec- 
tion criteria: 

• Impact parameter, with respect to the origin, in the r0 plane \Sr-(f,\ < 4.0 cm 

• Impact parameter, with respect to the origin, in the z plane \Sz\ < 6.0 cm 

• Icos^l < 0.94 
. ^ < 1.0 

• At least 1 r — (f) track hit from the silicon vertex detector (VD) 

• Tracks must not have been flagged as originating from material interactions (via 
any of the PXPHOT codes: -124, -123, -75, -85, -84 or -74). 

Neutral particles are used that are flagged by the DELPHI mass code as being a 
photon, 7r°, K° or A° (codes 21, 47, 61 and 81 respectively). 

2.3 Event Jets 

Event jets are reconstructed via the routine LUCLUS All reconstructed particles, 
that pass the standard selection, are used with a transverse momentum cutoff value of 
djoi„=PARU(44)=5.0. 
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2.4 Event Hemispheres 

All reconstructed particles, that pass the standard selection, are used to calculate the 
event thrust axis via routine LUTHRU Event hemispheres are then defined by the 
plane perpendicular to the thrust axis. 

Each hemisphere has associated with it two axes: the thrust axis and the refer- 
ence axis used for the calculation of rapidity (see Section |2.5| below) for tracks in that 
hemisphere. The reference axis is a jet axis selected in the following way: 

• If the event is a two jet event, the reference axis for the hemisphere is the jet axis 
in that hemisphere. 

• If a hemisphere contains 2 or more jets (i.e. an event with 3 or more jets) : 

— if one of the jets is the highest energy jet in the event, that jet axis forms the 
reference axis. 

— if the highest energy jet is in the opposite hemisphere, the combined probability 
(Pjet) for the tracks in a jet to have originated from the event primary vertex 
is formed (via the AABTAG package algorithms [^). 

The jet that is most 'B-like', i.e. with the smallest probability, is then selected 
if Pjet < 0.05. 

— if no jet in the hemisphere satisfies the above criteria, the jet with the highest 
energy is selected to form the reference axis. 

In addition, internally to BSAURUS, the hemispheres are numbered as 1 or 2 with 
the convention that the jet forming the reference axis for hemisphere 1 is of higher energy 
than that for hemisphere 2. 

2.5 Track Rapidity 

The rapidity of a track is defined as follows, 

for E the track energy and Pl the longitudinal momentum component along the reference 
axis for the hemisphere. 

The BSAURUS rapidity algorithm also returns an initial estimate of the B-hadron 
candidate 4-vector. This is defined simply to be the sum of individual track momentum 
vectors in a hemisphere for tracks with rapidity greater than 1.6. The value of the cut 
provides a good rejection of tracks originating from the primary vertex while accepting 
the majority of tracks from the B-decay. 

2.6 Particle Identification 

Unless otherwise stated, particle identification in BSAURUS is based on the following 
definitions: 

• Charged kaons: the kaon network output, or kaon probability, from the MACRIB 
package H and stored in BSAURUS output BSPAR(IBP_MACK,IH). 



4 



• Protons: the proton network output, or proton probability, from the MACRIB 
package Hand stored in BSAURUS output BSPAR(IBP_MACP,IH). 

• Electrons: the electron network output, or electron probability, returned by the 
ELEPHANT routine ELNETID [|T|. 

• Muons: the MUCAL2 definitions of a very loose, loose, standard and tight muon. 



2.7 Performance of Networks 

The main tool used to quantify and compare BSAURUS network performance is the 
purity against efficiency plot. Given some definition of 'signal' and 'background', purity 
is defined as, 

, correctly classified signal events 

lUrtty = (1) 

all classified events ' 



and the efficiency defined as, 

correctly classified signal events 
all signal events 



(2) 



In this note, the efficiency is calculated after all selection cuts, e.g. from Section |2.2| and 
Section ^ have been applied. When considering network outputs that involve symmetric 
binary decisions (i.e. the TrackNet and all fiavour tagging networks), we use an alternative 
measure of efficiency termed / = 'fraction classified' and defined as, 

^ all classified tracks /hemispheres 
all classifiable tracks /hemispheres 

For this case, in order to make the purity-efficiency plot, / is varied by making successive 
cuts in a band around the equal probability point for tagging 'signal' or 'background'. 



3 General BSAURUS Quantities 
3.1 Secondary Vertex Finding 

Per-hemisphere, an attempt is made to fit a secondary vertex to tracks with rapidity 
> 1.6 that pass the standard track selection criteria of Section P?^. 



From this class of track, additional criteria are applied with the aim of selecting tracks 
for the vertex fitting stage that are likely to have originated from the decay chain of a 
weakly decaying B-hadron state. As Figure |l| illustrates, it is important to reject as far as 
possible tracks from the fragmentation in order to avoid large pulls in the vertex position 
toward the primary vertex. There must be a minimum of two tracks for the vertex fit to 
be attempted and the selection process is as follows: 

[1] The highest energy muon or electron candidate is selected if Eiepton > 3 GeV. 

[2] The crossing point of each track with the initial B-candidate direction vector (see 
Section p75|) is found and the distance, L, between primary vertex [] and crossing 



^Unless otherwise stated, the primary vertex definition in BSAURUS is the event primary vertex from 
the AABTAG package. 
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Figure 1: The pull on the decay length residual L{recon) — L{truth) introduced if frag- 
mentation tracks are allowed in the secondary vertex fit definition. The open histogram 
shows the case where no tracks from fragmentation are in the fit and the hatched his- 
togram the case where fragmentation tracks are included. 
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point is calculated. Tracks are selected if L/a^ > 2 or, for cases where L/ai < 2, 
tracks must also satisfy y > 2.5 and L > 0.1cm. 



If at this point the number of tracks selected for the fit is less than 2, further attempts 
are made to add tracks to the fitting procedure. When this happens in bb events, it is 
likely that the true decay length of the B-hadron was small. 
If only one track was added from stages 1 and 2: 

[3] Add tracks with y > 3. 

[4] If still only one track selected, add the track of highest rapidity from the remaining 
track list of the hemisphere. 

If no track was selected from stages 1) and 2): 

[5] A search is made for the best kaon candidate in the hemisphere based on the kaon 
probability defined in Section If this kaon candidate has ?/ > 2 it is added to that 



track in the remaining track list of the hemisphere that has the highest rapidity. This 
track pair alone will then form the starting track list for the vertexing procedure. 

[6] Finally if no kaon candidate exists in the hemisphere, the two tracks of highest 
rapidity are selected. 

Using the track list supplied by this track selection process, a secondary vertex fit is 
performed in 3-dimensions (using routine DAPLCON from the ELEPHANT package ^) 
constrained to the direction of the B-candidate momentum vector. The event primary 
vertex is used as a starting point and if the fit did not converge^, the track making the 
largest contribution is stripped away in an iterative procedure, and the fit repeated. 

In addition to returning the secondary vertex position, the procedure also fits the 
primary vertex position and updates the B-candidate direction according to the vector 
joining the primary and secondary vertex points. This information is used in forming 
some of the track net inputs as described in the next section. 

Once a convergent fit has been attained, the final stage of the secondary vertex fitting 
procedure involves an attempt to add into the fit tracks that failed the initial track 
selection criteria but nevertheless are consistent with originating from the vertex. These 
tracks are identified on the basis of an intermediate version of the TrackNet. This is 
a neural network output that discriminates between tracks originating from the primary 
vertex and those likely to have come from a secondary vertex and is described in Section |^. 
A version of the TrackNet is constructed, specifically for the purpose of use in this final 
stage of vertex fitting, based on secondary vertexing information available before this final 
stage has run. In general, secondary vertex tracks have TrackNet ~ 1 whereas tracks from 
the primary vertex have TrackNet values close to zero. The track of largest TrackNet in 
the hemisphere is added to the existing track list and retained if the resulting fit converges. 
This process continues iteratively for all such tracks with TrackNet > 0.5. 

^Here, non-convergence means the fit took more than 20 iterations. A further iteration is deemed 
necessary if the is above 4 standard deviations during the first fO iterations or above 3 standard 
deviations during the next f iterations. 
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3.2 B-Candidate Four- Vectors Available from BSAURUS 



BSAURUS provides currently four different definitions of the B-candidate four- vector in 
a liemispfiere: 



The initial estimate from the rapidity algorithm outlined in Section |2.5| and stored 
in the four BSAURUS output array elements starting at BSHEM(IBH_BY, IH) for 
hemisphere number IH. 

The three momentum vector is then transformed to lie along the direction join- 



ing primary to secondary vertex as mentioned in Section The energy com- 

ponent is the rapidity energy and this vector is stored in location starting at 
BSHEM(IBH_BF,IH). 

In location starting at BSHEM(IBH_BT, IH), the energy component is then corrected 



to account for neutral energy losses (see Section [37^ ). The momentum vector derives 
from a TrackNet weighted momentum sum for the case of 2-jet events and is the 
rapidity vector for > 2-jet events. 

The four- vector starting at location BSHEM(IBH_BTT, IH) is, in principle, the same 
as the BSHEM(IBH_BT, IH) vector but contains a more detailed energy correction. 



described in detail in Section |3.3| . In addition, the momentum vector derives from 



the rapidity algorithm for > 2-jet events that also pass Ehem/Ebeam > 0.6 and 
from the TrackNet weighted sum otherwise. 



3.3 B-Candidate Energy Correction 

We describe here the energy correction procedure that went into making vector 
BSHEM(IBH_BTT,IH). (Note that the BSHEM(IBH_BT, IH) vector contains a very similar 
correction which, although still used internally inside BSAURUS, is older and less opti- 
mal.) 

Separate correction functions are derived for 2-jet and > 2-jet events. Within these two 
classes, because of the explicit energy ordering of the hemispheres, we treat hemispheres 
1 and 2 completely separately. To be used in the correction procedure, hemispheres must 
pass the following cuts: 

• The initial reconstructed B-candidate energy E^aw is 20GeV or more. 

• The corresponding initial reconstructed B-candidate mass rriraw lies within two stan- 
dard deviations of the total data sample median value. 

• The ratio, Xh, of the hemisphere energy Ehem to beam energy Eheam lies in the range 
0.6 <Xh< 1.1. 

In addition, for the case of hemispheres with larger missing energy, i.e. Xh < 0.6, a further 
correction is derived separately for hemispheres 1 and 2. 

The starting point is the initial estimates of the B energy and mass, E^aw and Tfi^a^. 
Following studies, we choose these estimates to be from the rapidity algorithm for events 
with > 2-jets and where Ehem/Ebeam > 0.6, and to be derived from the sum of 'B- 
weighted' four-vectors otherwise. This involves weighting (via a sigmoid threshold func- 
tion) the momentum and energy components of charged tracks by the TrackNet output 
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(see section ^ and those of neutrals by their rapidity. In this way the effects of tracks from 
the B-decay are enhanced and those of tracks from the primary vertex are suppressed. 

The correction procedure is motivated by the observation in Monte Carlo of a correla- 
tion between the energy residuals AE = E'^'^ — E^£'^^^ and mraw (which is approximately 
linear in rriraw), and a further correlation between AE and Xh resulting from neutral en- 
ergy losses and inefficiencies. The correction proceeds in the following way: the data are 
divided into several samples according to the measured ratio Xh and for each of these 
classes the B-energy residual AE is plotted as function of rriraw The median values of 
AE in each bin of rriraw are calculated and their rriraw dependence fitted by a third order 
polynomial 

AE{mraw; Xh) = a + h{mraw - (rriraw)) + C{mraw - {rrirawW + d{mraw " ("^ra^))^ 

The four parameters a, 6, c, in each Xh class are then plotted as a function of Xh and their 
dependence fitted with third and second-order polynomials. Thus one obtains a smooth 
correction function describing the mean dependence on rriraw and on the hemisphere 
energy as determined from the Monte Carlo. Finally, a small bias correction is applied 
for the remaining mean energy residual as a function of the corrected energy. 

Figure ^ shows typical resolutions on the reconstructed B-direction and energy attain- 
able from the four-vector output at location BSHEM(IBH_BTT, IHEM) . 

3.4 Decay Length and Decay Time 

The B-hadron decay length estimate in the r — plane, Lr^f), is defined as the (positive) 
distance between the primary and secondary vertex positions reconstructed from the 
secondary vertex search described in Section The three-dimensional decay length 
is constructed as L = Lr^/ sin(6'), where the B-hadron candidate momentum vector, pb 
detailed in Section |373| , defines the direction. The resulting proper decay time estimate is 
given by r = Lrris/pBC, where the B rest mass is taken to be 5.2789 GeV/c^. 

3.5 Track and Hemisphere Quahty 

When presenting neural networks with track and hemisphere-based discriminating vari- 
ables, it is advantageous at the network learning stage to provide the network with ad- 
ditional information relating to the potential quality of the input information. For this 
reason BSAURUS defines both track and hemisphere quality words that are used exten- 
sively in many network definitions. 

The track quality word stored in BSPAR(IBP_QUAL, IT) for track IT is defined as fol- 
lows: 

• 0, means the track is good. 

• 1, means the track contains some ambiguous track hit information fiagged by the 
DELPHI track hit ambiguity word. 

• 10, means that the track has been fiagged as originating from an interaction, defined 
by PXPHOT code=-120 returned by routine PXGECO. 

• 20, means that the track did not pass the internal track quality criteria of the 
AABTAG package. 
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Figure 2: For 2-jet events, (a) shows the residual between reconstructed and Monte Carlo 
truth for the B-direction and (b) the resolution of the reconstructed corrected-energy for 



hemispheres inside the selection cuts described in Section 3.3 
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The hemisphere quahty word, stored in, BSHEM(IBH_QUAL, IH) for hemisphere IH, has 
the following definition: 

• X 1, where N is the number of tracks in the hemisphere rejected by the standard 
track cuts (described in Section p.2|) 

• + X 100, where N is the number of tracks in the hemisphere flagged as originating 
from an interaction in the same way as for the track quality word 

• + X 1000, where N is the number of tracks in the hemisphere flagged as being 
ambiguous in the same way as for the track quality word 

• + X 100000, where A^ is the number of tracks in the hemisphere flagged as failing 
the quality cuts of the AABTAG package. 

For use in network training, the hemisphere quality is better used in a real, continuous 
form rather than as an integer. Stored in BSHEM(IBH_QUAL2, IH) therefore, is the trans- 
formation of the integer word into a continuous variable ranging from O.0(good quality) 
to 10.0(bad quality). 



4 Use of Neural Networks in BSAURUS 

Events and event hemispheres used to prepare the training samples for BSAURUS net- 
works must pass the following set of requirements: 

• 2-jet events with the jets back-to-back to better than 10°, 

• \cos6thrust\ < 0.65, 

• combined AABTAG event b-tag > 0.02, 

• both gas and liquid parts of the RICH must be functioning, 

• defining Eh em as the total energy contained in charged and neutral particles in the 
hemisphere, then a hemisphere must satisfy 0.6 < Ehem/Ebeam < 1-1, 

• a hemisphere must have an initial B-energy estimate from the rapidity algorithm 
> 20GeV, 

• only hemispheres with a fully converged secondary vertex fit are considered. 

BSAURUS uses the JETNET^] neural network package run in a feed-forward mode 
nearly always with three layers and a sigmoid function as the input-output (or transfer) 
function. The hidden layer typically contains one more node than the number of inputs. 

The network is trained using two statistically independent samples; the training sample 
and the performance sample. Usually each sample consists of an equal number of 'signal' 
and 'background' hemispheres/tracks. The training sample is used to teach the network 
to separate signal from background, but it is the network output error^ based on the 

■^Defined as \/^i{NOi — TARi)'^, where index i labels the network output node, NOi is the network 
output value at output node i and TARi is the target value for output node i. 
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performance sample that determines when the end-point of the training procedure has 
been reached. 

Variables used as inputs to the networks are selected based on (a)the discriminating 
power between signal and background and (b)the agreement between simulation and data. 
Variables are often also transformed and/or their range truncated so that the network is 
presented with variables that show, ideally, a linear, continuous variation in the ratio of 
'signal' to 'background' across the full range. 



5 The B-Track Probability Network (TrackNet) 

The aim of the B-track probability network is to discriminate between the class of tracks, 
in a bb-event hemisphere, originating from the weak decay of a B-hadron from all other 
tracks. 



Only tracks passing the standard quality cuts of Section 2.2 are considered. The 



discriminating variables that form the inputs to the network per track are: 

• The track total momentum. 

• The track momentum in the B-candidate rest frame. 

• The helicity angle of the track defined as the angle between the track vector in the 
B-candidate rest frame and the B-candidate momentum vector in the lab frame. 

• A flag to identify whether the track was in the secondary vertex fit or not. 

• The probability that the track originates from the fitted primary vertex (AABTAG 
algorithm) . 

• The probability that the track originates from the secondary vertex (AABTAG 
algorithm) . 

• The probability that the track originates from the fitted primary vertex (BSAURUS 
algorithm) . 

• The probability that the track originates from the secondary vertex (from the 
BSAURUS algorithm). 

• The track rapidity. 

In addition, input variables that gave no inherent discriminating power were included to 
inform the network of the potential quality of the other input variables: 

• The decay length or distance between the primary and secondary vertex in the r — 
plane. This quantity is then scaled by the reconstructed error to form a decay length 
significance. 



The track quality word from Section 3.5 



12 




12000 - 




Figure 3: Distributions of the input variables to the TrackNet, compared to data. Tracks 
are selected from events that are contained in the barrel, with 2-jets and are about 80% 
pure in T? — > bb events. The key to the symbols is given in the top left plot. Here, 'signal' 
refers to tracks originating from the B-hadron decay chain, 'background' are tracks from 
fragmentation or excited B-hadron decay, and u,d,s,c background are any tracks in non-b 
decays of the Z°. 
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Distributions of all discriminating variables are shown in Figure |^. Two points should 
be noted. Firstly, to first approximation the track probabilities to originate from the 
primary or secondary vertex as calculated by the BSAURUS or the AABTAG packages, 
are calculations of the same quantity. There are differences however in the way track 
errors are handled in the two cases and, in addition, the AABTAG package returns a 
probability tuned to the resolution seen in the data. For these reasons, the correlation in 
these variables between the two methods is not 100% and there is extra information to be 
gained by including all definitions as input variables to the network. Secondly, the decay 
length variable (which is the same for all tracks in a given hemisphere) and the track 
quality word do not a priori distinquish between B-decay products and the background. 
These variables do however provide the network with information about whether the 
hemisphere has a well separated secondary vertex and whether the information associated 
with a particular track is reliable or not. This additional information has been found to 
be useful during the network training procedure. 

The network therefore consists of 11 input nodes, one for each of the variables listed 
above, 12 nodes in the hidden layer, and a single output node with target value 1 for 
signal tracks and for background. The training and performance samples are drawn 
only from Z° bb Monte Carlo where signal tracks are defined as coming from any point 
along the weak B-decay chain, and background as coming from any source located at 
the primary vertex e.g. fragmentation and excited B-hadrons. The network training 
was performed using 30K signal and 30K background tracks in both the training and 
performance samples. The network output and performance for tagging tracks from B- 
decay are shown in Figures ^ and |^. 

6 B-Species Identification(SHBN) 

BSAURUS returns in the output variables, BSHEM(IBH_PRBBP,IH), 
BSHEM(IBH_PRBBO,IH), BSHEM(IBH_PRBBS , IH) and BSHEM(IBH_PRBLB , IH) the probabil- 
ity for the hemisphere IH to contain a B"*", B", Bs or B-baryon respectively. Q The approach 
employed is to supply all discriminating variables to a neural network with four output 
nodes, one for each B-hadron type. 

The resulting Same Hemisphere B-species enrichment Network or SHBN, consists of 
15 input nodes(described below) and 17 nodes in the single hidden layer. Each output 
node delivers a probability for the hypotheses it is trained on i.e. the first supplies the 
probability for a Bg meson to be produced in the hemisphere, the second for a B° meson, 
the third for charged B"*" mesons and the fourth for all species of B-baryons. The following 
input variables are used: 

• Using the TrackNet value as a probability (Pb) for each track to originate from the 
B-hadron decay vertex rather than from the primary vertex, the weighted vertex 
charge is formed 

tracks 

\Q,\= Pb{^-Q{^ (4) 

i 

that distinguishes between charged and neutral B-hadrons and is illustrated in Fig- 
ure 1^. 

''The charge conjugate states are also implied. 
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Figure 4: The TrackNet output distribution for Monte Carlo compared to the data. 
Tracks are selected from events that are contained in the barrel, with 2-jets and are 
about 80% pure in Z° bb events. The shaded distributions are defined as in Figure |^. 
Here, 'Signal' refers to tracks originating from the B-hadron decay chain, 'background' 
are tracks from fragmentation or excited B-hadron decay, and u,d,s,c background are any 
tracks in non-b decays of the Z°. 
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Figure 5: The TrackNet performance. Tracks are selected from events that are contained 
in the barrel, with 2-jets and are about 80% pure in Z° bb events. The plot is based 
on equal numbers of 'signal' tracks that originate from the B-hadron decay chain, and 
'background' tracks originating from fragmentation or excited B-hadron decay. 
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Figure 6: The vertex charge variable comparing data(points) to simulation(histogram) 
and showing separately the contributions from neutral, positively and negatively charged 
B-hadrons. Selection cuts and weights have been applied as detailed in Appendix A. 
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• The binomial error on the vertex charge, defined as 

= y: ^PBm-PB{i)) (5) 

i 

gives a measure of the rehabihty of the vertex charge information. 

• The number of charged pions in the hemisphere. This is most powerful for the case 
of B-baryons and Bg mesons, which have a higher content of non-pion particles, i.e. 
neutrons, protons and kaons in comparison to other B-species. 

• The total energy deposited in the hemisphere by charged and neutral particles 
scaled by the LEP beam energy. This is sensitive to the presence of B-baryons 
and Bg mesons, due to the fact that associated neutrons and K^i^ are often not re- 
constructed in the detector, with the consequence that the total hemisphere energy 
tends to be smaller compared to B'' or B+ mesons. 

• Bs mesons are normally produced with a charged kaon as leading fragmentation 
particle with a further kaon emerging from the weak decay (the same applies to 
the associated production of protons with B-baryons). Exploiting this fact, input 
variables are constructed giving the likelihood for the presence of a leading frag- 
mentation kaon/proton and a kaon/proton weak decay product, for each of the two 
hypotheses; Bg or B-baryon. 

Utilising the Monte Carlo truth information, normalised track rapidity distributions 
were parameterised with Gaussians separately for the case of leading fragmentation 
tracks and B-decay products. These were then used to form a weight per recon- 
structed track which was summed over to give a hemisphere level variable. Cuts were 
applied for kaon or proton identification (see Section |2.6|) and in order to separate 
the fragmentation tracks from the B-decay tracks via the TrackNet output. 

• The leading fragmentation track also can often be neutral e.g. production of a 
Kg associated with Bg meson production or A° associated with B-baryon production. 
An input was therefore constructed based on the presence of a reconstructed K° or 
A*^ in the hemisphere in the same way as for charged kaons and protons described 
above. In this case, the rapidity distribution was formed from Monte Carlo truth 
information for neutral particles originating from the primary vertex. The resulting 
weights were then summed over for reconstructed neutral particles with DELPHI 
mass codes 61 or 81 identifying K° and A° candidates. 

• The probability for the leading fragmentation particle to be a charged kaon. 
Specifically, the maximum kaon net output from the three tracks with highest 
rapidity originating from the primary vertex (via the condition TrackNet< 0.5) 
was used. The three kaon probabilities can be found in BSAURUS outputs, 

BSHEM ( IBH.KRANKl , IH) , BSHEM (IBH_KRANK2 , IH) , BSHEM ( IBH_KRANK3 , IH) . 

^The term leading refers to the neighbouring hadron to the B-hadron that emerges from the fragmen- 
tation chain. This particle is identified in practice from the fact that it often has the largest rapidity with 
respect to the B-hadron flight direction. 
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• A charge correlation between the leading fragmentation particle charges provided a 
flag for the presence of B"*" mesons. Specifically, the rapidity-weighted track charge 
sum over all tracks in the hemisphere is formed, and this is then scaled by the 
measured vertex charge. 

In addition, input variables that gave no inherent separation power between different B 
species were included to inform the network of the potential quality of the other input 
variables: 

• the invariant mass of the reconstructed vertex, 



the hemisphere quality factor described in Section 



• the energy of the B-hadron, to provide information on how hard the fragmentation 
was and therefore inform the network of how the available energy is expected to be 
shared between B hadron and fragmentation products. 

The output of the SHBN in both simulation and data is shown in Figure ^ and the 
resulting performance of the SHBN is shown in Figure |^. The performance can be 
further improved by the combination of charge correlation information from the opposite 
hemisphere to form the Both Hemispheres B-species enrichment Network or BHBN. This 
network is described in Section p.3.1| . 



7 B-D Separation(BDnet) 



In selecting tracks for inclusion in the B-secondary vertex fit described in Section |3.1| , 
there is inevitably some background from tracks that originate not from the B-decay 
vertex directly, but from the subsequent D-cascade decay point. The effect of including 
such tracks is illustrated in Figure § which shows that, on average, the B-decay length 
estimate reconstructed by BSAURUS lies somewhere between the true B decay point and 
that of the cascade D. 

The BDnet is a neural network designed to discriminate between tracks originating 
from the weakly decaying B-hadron and those from the subsequent cascade D-meson de- 
cay. This information is clearly crucial to any algorithm wishing to eliminate the positive 
decay length bias just described, and is also used extensively throughout BSAURUS, 



particularly in the construction of B- flavour tags (e.g. see Section 
The following discriminating variables are used: 

• The angle between the track vector and an estimate of the B flight direction (taken 
to be BSHEM(IBH_BTT,IH)). 

• The probability that the track originates from the fitted primary vertex (AABTAG 
algorithm) . 

• The probability that the track originates from the fitted secondary vertex (AABTAG 
algorithm) . 

• The momentum and angle of the track vector in the B rest frame. 

• The TrackNet output. 
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Figure 7: Output of the SHBN network for the B"'",B°,Bs and B-baryon hypotheses in 
the simulation(histogram) compared to the data(points). The component histograms in 
each case show the distribution for the hypothesis being considered, i.e. the 'signal' (open 
histogram), compared to the distribution for everything else, i.e. the 'background' (shaded 
histogram). Selection cuts and weights have been applied as detailed in Appendix A. 
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Figure 8: Performance of the SHBN for enriching samples in the various B-hadron types. 
The plot is based on bb Monte Carlo making successive cuts in the network output and 
where 'signal' and 'background' are defined as in Figure ^. Selection cuts and weights 
have been applied as detailed in Appendix A. 
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Figure 9: Comparing {L{recon) — L{B)) to {L{recon) — L{D)) illustrating that the 
reconstructed decay length is often larger than the true B-decay length but less than the 
true D-cascade decay length. The origin of this pull to larger decay lengths is the inclusion 
of one or more tracks from the D-vertex. 
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The kaon network output, described in Section O. 



The lepton identification, see Section ^ 



The remaining variables carry no imphcit discriminating power but are included as gauges 
of the quality of the other variables: 

• The track quality word from Section T^. 

• The hemisphere quality word from Section |3.5| . 

• The hemisphere decay length significance, L/ai- 

• The hemisphere secondary vertex mass. 

• The hemisphere rapidity gap between the track of highest rapidity below a TrackNet 
cut at 0.5 and that of smallest rapidity above the cut at 0.5. 

The network architecture consisted of 14 input nodes with 15 nodes in a hidden layer and 
with a single output node trained on equal numbers of tracks directly from B-decay as 
'signal' or from the following D-decay as 'background'. Figure [1^ shows the output of the 
network for the two classes of track used in training and Figure ^ plots the performance 
of the tag. 



8 Optimal Fragmentation, Decay Flavour and Decay 
Type Tagging 

In an attempt to use the information available optimally, the BSAURUS approach to 
tagging the fiavour, i.e. charge of the b-quark, works by first constructing a probability 
per track and then combining these to give a probability at the hemisphere level. 

Separate neural networks are trained to tag the underlying quark charge for the cases 
when there was a B+,B°,Bs or B-baryon produced in the hemisphere. In addition two 
sets of such networks are produced, one trained only on tracks originating from the frag- 
mentation process and the other trained only on tracks originating from weak B-hadron 
decay. 



8.1 The Basis of Optimal Flavour Tagging 

The basis of optimal fiavour tagging is to construct, at the track level, the conditional 
probability for the track to have the same charge as the b-quark in the B-hadron both 
at the moment of fragmentation (i.e. production) and at the moment of decay. In ad- 
dition, these probabilities are constructed separately for each of the B-hadron types; 
B+,BO,Bs and B-baryon. 

A neural network is used with a target output value of +1(-1) if the track charge is 
correlated(anticorrelated) with the b-quark charge. The discriminating input variables 
used are the following: 

• Particle identification variables (see Section P^ ): 
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Figure 10: Output of the BDnet for tracks in simulation(histogram) and in data(points), 
for tracks with TrackNet> 0.5. The component histograms show the two classes of track 
that the network was trained on, namely tracks originating from cascade D-decay('signar) 
and all other tracks which are mainly tracks from B-decay('background'). Tracks are 
selected from events that are contained in the barrel, with 2-jets and are about 80% pure 
in Z° bb events. 
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Figure 11: The purity against efficiency performance of tiie BDnct. Tlie plot is based 
on equal numbers of 'signal' tracks, i.e. B-D cascade tracks, and 'background' tracks, i.e. 
all other tracks with TrackNet> 0.5. Tracks are selected from events that are contained 
in the barrel, with 2-jets and are about 80% pure in Z° bb events. 
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— kaon network output, 

— proton network output, 

— electron network output, 

— muon classification code. 

B-D vertex separation variables: 

— BDnet output, 

— ^'^"^j^™"^ ; where BD is the BDnet value, BD{min) is the minimum BDnet 
value for all tracks in the hemisphere above a TrackNet value of 0.5, and 
ABD is the difference between BD{max) and BD{min), 

— the momentum of the track in the B-candidate centre of mass frame. 
Track-level quality variables: 

— the helicity angle of the track vector in the B-candidate centre of mass frame. 



— the track quality as defined in Section ^3 

— the TrackNet output, 

— the track energy. 

• Hemisphere-level quality variables: 

— the hemisphere rapidity gap between the track of highest rapidity below a 
TrackNet cut at 0.5 and that of smallest rapidity above the cut at 0.5, 

— the hemisphere quality as defined in Section 

— the number of tracks in the hemisphere passing the standard cuts (Section |2.2|) 
and with a TrackNet value greater than 0.5, 

— ABD as defined above, 

— the secondary vertex mass, 

— the secondary vertex fit probability, 

— the ratio of the B-energy, defined by the BSHEM(IBH_BTT, IH) , to the LEP beam 
energy, 

— the error on the vertex charge measurement. 

In total, the track decay flavour network used all input variables described above in a 
network containing two hidden layers with 20 nodes in the first and 3 in the second. 
The track production flavour network used the same input variables without the lepton 
identification and B-D vertex separation variables, with 19 nodes in the first hidden layer 
and 3 in the second. 

The resulting track charge correlation conditional probabilities are plotted in Fig- 
ure O for the case of fragmentation flavour and in Figure O for the case of de- 



cay flavour. Figures and ^ show the corresponding tag performance attained 
by the fragmentation and decay flavour probabilities respectively. These outputs 
can be found in BSAURUS variables; BSPAR(IBP_FFLBS, IT),BSPAR(IBP_FFLBO , IT), 
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Figure 12: The track-level conditional fragmentation flavour probability for all B-liadron 
types, comparing data(points) and simulation(histogram). The plots with larger normal- 
isation correspond to the normal mixture of B-hadron types whilst the distibutions with 
smaller normalisation are for samples enhanced in that B-type. The enhancement cuts 
were chosen to give roughly a 10% selection efficiency of signal, i.e. BHBN greater than: 
0.8 for B^, 0.6 for B*^, 0.2 for Bg and 0.33 for B-baryons. All tracks have TrackNet less 
than 0.5 and the events used were 2-jet events only, contained in the barrel. Weights were 
applied as in Appendix A 
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Figure 13: The track-level conditional decay flavour probability for all B-liadron types, 
comparing data(points) and simulation(histogram). The plots with larger normalisation 
correspond to the normal mixture of B-hadron types whilst the distibutions with smaller 
normalisation are for samples enhanced in that B-type. The enhancement cuts were 
chosen to give roughly a 10% selection efficiency of signal, i.e. BHBN greater than: 0.8 
for B+, 0.6 for B°, 0.2 for Bg and 0.33 for B-baryons. All tracks have TrackNet greater 
than 0.5 and the events used were 2-jet events only, contained in the barrel. Weights were 
applied as in Appendix A 
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BSPAR(IBP_FFLBP,IT) and BSPAR(IBP_FFLLB, IT) for the fragmentation 
flavour conditional probabilities and in BSPAR(IBP_DFLBS, IT) ,BSPAR(IBP_DFLBO , IT) , 
BSPAR(IBP_DFLBP, IT) and BSPAR(IBP_DFLLB, IT) for the decay flavour conditional prob- 
abilities. 

Finally, to obtain a flavour tag at the hemisphere level, the conditional track prob- 
abilities described above, P{sameQ\iy where i = B"^,B°,Bs or B-baryon and j = 
fragmentation or decay, are combined as the likelihood ratio, 

F[hem.)\ = 2^ In — — — .Q{track) 

where Q{track) is the track charge. Which tracks to sum over depends on the hypothesis 
considered, i.e. for the fragmentation flavour it is over all tracks with TrackNet< 0.5 whilst 
for the decay flavour, tracks must satisfy TrackNet> 0.5. These hemisphere flavour tags 



are compared with the data in Figure |Tg for the fragmentation tag, and in Figure |T7| for 
the decay tag. 



The purity against efficiency performance of the tags are shown in Figure |T8| for 
the fragmentation tag and in Figure for the decay tag. The tags are located 
in the following BSAURUS variables; BSHEM(IBH_FFLBS , IH) , BSHEM(IBH_FFLBO , IH) , 
BSHEM(IBH_FFLBP,IH) and BSHEM(IBH_FFLLB, IH) for the production flavour out- 
puts and in BSHEM(IBH_DFLBS,IH), BSHEM(IBH_DFLBO,IH), BSHEM(IBH_DFLBP, IH) and 
BSHEM-[IBH_DFLLB,IH) for the decay flavour outputs. 

8.2 Applications to Production Flavour Tagging 

8.2.1 The Same Hemisphere Production flavour Network (SHPN) 

The hemisphere flavour tags described in Section form the major input to a further 
dedicated network, the Same Hemisphere Production flavour Network or SHPN. This 
network attempts to find the optimal combination of fragmentation and decay flavour, 
for each B-species hypothesis, in order to tag the B-hadron quark flavour at production 
time. The network is constructed independently of any information from the opposite 
hemisphere, primarily so that the output can be used for a ApB{bb) measurement |TD[ 
incorporating double-hemisphere flavour tag methods. 

The 9 input variables used are (ignoring some details of variable transformation and 
re-scaling) : 

• FQiem^^^l"-^- ■ Psame{Bs) 

. (F(/iem.)gf ■ - F{hem.)lT-) ■ Psan.e{B^) 

• (^F{hem.)j^^;y^ - F(/iem.)f„7^^-) ■ Psame{hary.) 
F{hem.)'^r- ■ (l - 2sin2(0.237 ■ r)) - F(/iem.) Jji"^-) ■ P,ame(5°), where r is the re- 



constructed B-lifetime. Note this construction takes account of the B*^ oscillation 
frequency found in the simulation. This is not possible for the case of Bg where the 
oscillations are so fast that we have essentially a 50-50 mix of Bg and Bg. 
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Figure 14: The track-level production flavour tag performance. All tracks have TrackNet 
less than 0.5 and the events used were 2-jet events only, contained in the barrel. 
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Figure 15: The track- level decay flavour tag performance. All tracks have TrackNet 
greater than 0.5 and the events used were 2-jet events only, contained in the barrel. 
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Figure 16: The hemisphere-level fragmentation flavour tag for all B-hadron types, com- 
paring data(points) and simulation(histogram). The plots with larger normalisation cor- 
respond to the normal mixture of B-hadron types whilst the distibutions with smaller 
normalisation are for samples enhanced in that B-type. The enhancement cuts were cho- 
sen to give roughly a 10% selection efficiency of signal, i.e. BHBN greater than: 0.8 for 
B+, 0.6 for B'^, 0.2 for Bg and 0.33 for B-baryons. Selection cuts and weights have been 
applied as detailed in Appendix A. 
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Figure 17: The hemisphere-level decay flavour tag for all B-hadron types, comparing 
data(points) and simulation (histogram). The plots with larger normalisation correspond 
to the normal mixture of B-hadron types whilst the distibutions with smaller normalisa- 
tion are for samples enhanced in that B-type. The enhancement cuts were chosen to give 
roughly a 10% selection efficiency of signal, i.e. BHBN greater than: 0.8 for B+, 0.6 for 
B°, 0.2 for Bg and 0.33 for B-baryons. Selection cuts and weights have been applied as 
detailed in Appendix A. 
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Figure 18: The hemisphere-level fragmentation flavour tag performance for all B-hadron 
types. Selection cuts and weights have been applied as detailed in Appendix A. 
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Figure 19: The hemisphere- level decay flavour tag performance for all B-hadron types. 
Selection cuts and weights have been applied as detailed in Appendix A. 
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• The jet charge defined as, 

where the sum is over all tracks and is the longitudinal momentum component 
to the thrust axis. The optimal choice of the free parameter k depends on the type 
of B-hadron under consideration. For this application we choose a range of values 
K = 0.3, 0.6, oo where the last value corresponds to taking the charge of the stiffest 
track in the hemisphere. 

• Vertex charge defined as in Section ^ equation ^. 

• Vertex charge significance. 

Note that the Psame{B~^ , B^, Bs, hary.) factors represent the outputs of the B-species tag- 
ging network, SHBN, described in Section ^. 

The treatment of quality variables is slightly different for this network compared to 
other BSAURUS networks. In this case, in order to ensure that the output is inherently 
symmetric with respect to opposite charges, the quality variables are used to weight the 
turn-on gradient (or 'temperature') of the sigmoid function used as the network node 
transfer function. This restriction is especially important for analyses such as measure- 
ments of ApBibb) which to first order implicitly assume that the charge tag used is sym- 
metric with respect to quark and anti-quark. The quality variables used were: 



The hemisphere quality word BSHEM(IBH_QUAL2, IH) described in Section 3.5 



• The hemisphere rapidity gap between the track of highest rapidity below a TrackNet 
cut at 0.5 and that of smallest rapidity above the cut at 0.5. 

• The error on the vertex charge measurement. 

• The ratio of the B-energy, defined by the BSHEM(IBH_BTT, IH) , to the LEP beam 
energy. 

The SHPN output is stored in BSAURUS variable BSHEM(ISH_FLAV45 , IH) 
and is compared to the data in Figure The performance of the tag can be seen in 
Figure ^(full circles) . 



8.2.2 The Both Hemispheres Production flavour Network (BHPN) 

It should be noted that an improved performance for production fiavour tagging can be 
obtained by including the information from the opposite side hemisphere. The Both 
Hemispheres Production fiavour Network or BHPN is defined as the ratio, 

SHPN(same) ■ (1 - SHPN(oj9j9)) 
SHPN(same) ■ (1 - SHPN(opp)) + (1 - SHPN(same)) ■ SHPN(opp) ^ ^ 

where the labels same and opp refer to the production tag in this hemisphere and in 
the opposite hemisphere respectively. Note that if SHPN tags e.g. the charge of the 
quark, then (1 — SHPN) tags the charge of the anti-quark. The resulting performance 
of the BHPN is also shown in Figure |21|(open circles) and the output can be found in 
BSAURUS output BSHEM(IBH_FLBOTH,IH). 
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Figure 20: Output of the production flavour net in simulation compared to data. Over- 
laid are plots from the simulation showing the separation attained for the case where a 
produced B-hadron contains a b or b quark. The black region shows the distibution for 
the light and charm-quark background. Selection cuts and weights have been apphed as 
detailed in Appendix A. 
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Figure 21: Perfomance of the production and decay flavour networks. The comparison is 
made to using only the vertex charge or jet charge variable (with k = 0.6) as a production 
flavour tag. Selection cuts and weights have been applied as detailed in Appendix A. 



38 



8.3 Applications to Decay Type and Decay Flavour Tagging 

8.3.1 Both Hemispheres B-species enrichment Network (BHBN) 

The performance of the SHBN (Section ^ can be further improved by use of the hemi- 
sphere flavour tags developed in Section ^ and by the introduction of charge correlation 
information from the opposite hemisphere to form the Both Hemispheres B-species en- 
richment Network or BHBN. 

In the same way as for the SHBN, the BHBN is a network with four output target 
nodes, one for each B-type, 16 input nodes and two hidden layers consisting of 20 and 
10 nodes respectively. The input variables used, ignoring some details of variable scaling 
and transformation, are as follows: We define the variable, FLqpp, which is proportional 
to the production flavour tag SHPN(see Section p.2| ) in the opposite hemisphere, but 
scaled by a function in 1 — cos(6') in an attempt to explicitly account for the b-quark 
forward-backward asymmetry, 

FLopp = In {SHPN{opp)/{l - SHPN{opp))) - 
In ((0.5 - 0.12 ■ cos(^))/(0.5 + 0.12 ■ cos(^))) 

The discriminating variables used are: 

SHBN probability for B+, Psame{B+). 
SHBN probability for B^, Psan,e{B^)- 
SHBN probability for Bs, Psame{Bs). 
SHBN probability for B-baryon, Psameipary .) . 
F{hem.)^r- ■ {FLoPP - F{hem.f/:^-) 
FLopp ■ F(/iem.)^™'^ 



F{hem. 



Dec. 

Bs 


{FLqpp — 


F{hem.)pl°'^-) 




Dec. 
BO 


{FLqpp — 


F{hem.)lr-) 




Dec. 
bary. 


■ {FLqpp — 


F(/iem.)f„7/-) 




Dec. 

B+ 


{FLqpp — 


F{hem.)^+^') 


^ same \^ ) 


Dec. 

Bs 


{FLqpp — 


F{hem.)l:^^-) 


Psame {Bg) 


Dec. 
BO 


{FLqpp — 


F{hem.)lr-) 


same \^ ) 


Dec. 
bary. 


■ {FLqpp — 




■ Psame{bary 



F{hem. 

Where the F{hem.) terms are the hemisphere flavour likelihood ratios described in Sec- 



tion 8.1. The quality variables used are: 



The hemisphere quality word BSHEM(IBH_QUAL2, IH) described in Section ^3 
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The hemisphere rapidity gap between the track of highest rapidity below a TrackNet 
cut at 0.5 and that of smallest rapidity above the cut at 0.5. 

The ratio of the B-energy, defined by the BSHEM(IBH_BTT, IH) , to the LEP beam 
energy. 



The outputs of the BHBN are plotted in Figure ^ and can be found in BSAURUS 
variables, BSHEM(IBH_PRBAIBP, IH), BSHEM(IBH_PRBAIBO, IH), BSHEM(IBH_PRBAIBS, IH) 
and BSHEM(IBH_PRBAILB,IH). The resulting performance of the network is shown in fig- 

m. 



ure 



8.3.2 The Both Hemispheres Decay Flavour Network(BHDN) 

The Both Hemispheres Decay flavour Network or BHDN, utilises similar constructions 
for the discriminating variables as in the SHPN, the main difference being that opposite 
side hemisphere information is explicitly used in order to boost the performance. 

We define the variable, FLqpp, which is proportional to the production flavour tag 
in the opposite hemisphere, FL, but scaled by a function in 1 — cos(6') in an attempt to 
explicitly account for the b-quark forward-backward asymmetry, 

FLqpp = In {S H P N (opp) / {1 - SHPN{opp))) - In (1 - 0.4 ■ cos{e)) (8) 

The discriminating variables used are: 

. F(/iem.)gf • ■ PtothiBs) 

• (F(/iem.)gf ■ + FLqpp - F{hem.)lT-) ■ P,oth{B+) 
. (F(/iem.)f- + FLqpp - F{hem.)ZT) ' Pbothibary) 

• (F(/iem.)g?^- + {FLqpp - Fihem.)^^-) ■ (l - 2 sin2(0.237 * r))) -PbothiB"), where 
r is the reconstructed B-lifetime. 

Note that the Pboth{Bs, B^, bary, B^) factors represent the outputs of the B-species tagging 
network, BHBN which includes information from the opposite side hemisphere, and is 
described in Section |^. The quality variables used were: 

• Hemisphere quality (word BSHEM(IBH_QUAL2 , IH) described in Section ^3]) . 



The hemisphere rapidity gap between the track of highest rapidity below a TrackNet 
cut at 0.5 and that of smallest rapidity above the cut at 0.5. 



FLqpp. 

,2 



1 - 2sin^(0.237* r). 
BSHEM(IBH_FLAV45,IH), see Section 



The output of the decay flavour net is stored in BSAURUS word 
BSHEM(ISH_DFLAV45 , IH) and the performance of the tag is also shown in Figure ^(open 
squares) . 

BSAURUS routine, DECFLA, gives the user the option to recalculate the decay flavour 
network output based only on a list of tracks that he/she supplies. This can be useful in 
analyses that need to avoid potential biases in the flavour tag. 
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Figure 22: Output of the BHBN network for the B"'",B°,Bs and B-baryon hypotheses 
compared to the data (closed points). The two overlaid curves show the distribution 
for the hypothesis being considered i.e. the 'signal' (open histogram), compared to the 
distribution for everything else i.e. the 'background' (shaded histogram). Selection cuts 
and weights have been applied as detailed in Appendix A. 
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Figure 23: Performance of the BHBN for enriching samples in the various B-hadron 
types. The plot is based on bb Monte Carlo making successive cuts in the network output 



and where 'signal' and 'background' are defined as in Figure |2^. Selection cuts and weights 
have been applied as detailed in Appendix A. 
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Summary and Outlook 



The BSAURUS package reconstructs properties of B-hadron decays inside the jets of 
— > bb events in an intrinsically inclusive way in order to ensure high efficiency. Ex- 
ploiting wherever possible, physics knowledge of b-quark production, the hadronisation 
process and of the subsequent B decay, this information is combined to tag quantities of 
importance for B-physics analyses. In order to combine information in an optimal way, 
i.e. to ensure high purity, neural network techniques are used extensively. 

This note has described how this approach is applied to DELPHI data in order to 
form, most importantly, variables that tag the presence of a B"*", B°, Bs or B-baryon and 
variables that tag the charge or flavour of the b-quark in any of the B-hadron types at 
production and decay time separately. These tags are fundamental inputs to B-physics 
analyses such as B-spectroscopy0], B-species lifetime measurement s||], B° — Bj] and B° — 
B° oscillations 0, the forward-backward asymmetry in bb event s |llO[] , the measurement of 
B-hadron production fractions|rT| and measurements of B-hadron branching ratios . 
Analyses in all these areas have either already been completed or are currently underway in 
DELPHI using the BSAURUS framework. In addition, the list is by no means exhaustive 
and new analyses are currently being planned e.g. measuring the fragmentation function 
of B-species and studies of radially excited B-meson states. 

It should be noted that the package is under constant and on-going development. New 
features and improvements are still being added and this is likely to continue in the future. 
All news and new features relating to BSAURUS are documented at the web-site location: 
http : //pubxx . home . cern. ch/pubxx/tasks/bcteam/www/ inclusive/bsaurus .html . 
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Appendix A: Cuts and Weights Applied to BSAURUS 
performance Plots 

Where indicated, the following cuts have been applied to the event sample entering the 
distribution shown: 

• BSAURUS word, ISAURUS(IH), for hemisphere number IH must have value 4 or 
higher indicating that the secondary vertex fit in this hemisphere converged nor- 
mally. 

• Combined event b-tag> —0.5. 

• 2-jet events only. 

• \co's,{OTHmsT)\ < 0.75. 
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In addition, to account for known differences between data and Monte Carlo arising 
from both detector effects and physics modelhng, the simulation has been weighted (via 
BSAURUS routine SIMWAIT, stored in output BSHEM(IBH_WEIGHT, IH)) to account for: 

• The differences between data and Monte Carlo in the hemisphere quality vari- 
able BSHEM(IBH_QUAL, IH) as a function of the hemisphere track multiplicity, 
BSHEM(IBH_NPGOOD,IH). 

• Latest world average B-species production fractions. 

• Latest world average B-species lifetime measurements. 

• MARKIII D-topological branching fractions. 

• Any remaining difference between data and Monte Carlo in the B-energy estimate 
i.e. variable BSHEM(IBH_BTT+3,IH). 

Note that the current measurements of B-physics quantities used to form weights, are 
chosen to follow LEP Working Group recommendations. 



Figure ^ shows the effect, for example, on the production flavour network output of 
applying the weight returned by routine SIMWAIT. The much improved agreement between 
data and simulation after the weight is applied illustrates the need for such a correction 
if the network outputs are used in analyses to estimate e.g. absolute efficiencies. 



Appendix B: How to Access BSAURUS Outputs 

A call to routine BSAURUS should be made once per event from the user routine. A list 
of the BSAURUS outputs at track, hemisphere and event level is given in Appendix C. 

BSAURUS uses a modified version of the VECSUB package, called VECTOP, to fill 
the standard track information array P(NDIM,*). Of the utmost importance to users 
running the standard VECSUB package e.g. SKELANA users, is the track output word, 
BSPAR(IBP_IYB,IT). This gives the matching VECSUB track index for BSAURUS (ie 
VECTOP) track IT. Users failing to do this track matching will certainly obtain the 
wrong results from BSAURUS. 

The BSAURUS-cards file is released periodically to the standard DELPHI areas for 
library creation. A development version is to be found at: 
/af s/cern. ch/user/b/barker/public/delphi/bsaurus . car. 
Naturally, the results of using this version cannot be guaranteed. 



Appendix C: BSAURUS output array definitions 

Quantities calculated in BSAURUS along with other variables, e.g. truth Monte Carlo 
information, that are useful for analyses are available via COMMON /CSAURUS/ at the event, 
hemisphere and track level. 
1) The event output array: 

* Array Element Number of Description 

* values 
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Figure 24: The distributions show the ratio of data/simulation for the output of the pro- 
duction flavour network before(upper plot) and after(lower plot) applying the weighting 
procedure of Appendix A. 
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BSEVT(IBE_BTAG) 


1 


event B-tag 


BSEVT(IBE_BCONF) 


1 


event B-confidence (Bill Murray) 


BSEVT(IBE_BCTAG) 


1 


event B-tag with charm suppression 


BSEVT(IBE_BCOMB) 


1 


event aabtag combined-tag 


BSEVT(IBE_PV) 


3 


primary vertex 


BSEVT(IBE_PC) 


6 


primary covariance matrix 


BSEVT(IBE_NJET) 


1 


number of jets 


BSEVT(IBE_COST) 


1 


COS(theta) of event thrust axis 


RCJFVTfTRF RTCH") 


1 

J. 


sdclj-J-ksj. rtJ-on oucluu.o> o — KjCIo~xj.u. xj.u. uiixy 






,l=gas only, 0=no liq. or gas 


BSEVT(IBEMC_FLAV) 


1 


MC primary quark flavour (e.g. 12=b) 


BSEVT(IBEMC_PV) 


3 


MC primary vertex 


The hemisphere output array where IHEM is the hemisphere number (1 or 2): 


Array Element 


Number 


of Description 




values 





BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 



,HNUM,IHEM) 

BTAG.IHEM) 

BCONF.IHEM) 

BCTAGdHEM) 

BCOMB.IHEM) 

,B JET, IHEM) 

N JET, IHEM) 



BSHEM(IBH_BF,IHEM) 
BSHEM(IBH_BT,IHEM) 
BSHEM(IBH_BY,IHEM) 



BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 
BSHEMdBH. 



.BTT.IHEM) 
.SV,IHEM) 
.SCIHEM) 
.SIGPHI,IHEM) 
.SIGTHET,IHEM) 
.DEC, IHEM) 
.DECS , IHEM) 
.PROB,IHEM) 
.RET, IHEM) 
.YBFIT,IHEM) 
.NACC , IHEM) 
.Q JET, IHEM) 
.QJET2,IHEM) 
.Q JETS, IHEM) 
.QVER,IHEM) 
.QVERS , IHEM) 
.RAWQ , IHEM) 



1 
1 
1 
1 
1 
4 
1 

4 
4 
4 

4 
3 
6 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



hemisphere number 
hemisphere B-tag 

hemisphere B-confidence (Bill Murray) 
hemisphere B-tag with charm suppression 
hemisphere aabtag combined tag 
hemisphere B-candidate jet 4-vector 
number jets in hemisphere 

B 4-vector from fit 

B 4-vector from fit + energy correction 
B RAW 4-vector: from rapidity algorithm 
for 3-jet, net weighted for 2-jet 
'optimal^ B 4-vector 
secondary vertex coordinates 
secondary vertex covariance matrix 
Error on direction phi from SV fit 
Error on direction theta from SV fit 
B decay length 
B decay length error 
secondary vertex fit probability 
ISAURUS - fit return code 
YBFIT2 - secondary vertex fit return code 
no. of particles at sec. vertex 
jet charge, kappa=0.6 
jet charge, kappa=0.3 
jet charge, stiff est track in P_L 
Vertex charge NN 
Vertex charge NN error 
Vertex charge (sum trk. charge) 
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BSHEM( 


:iBH. 


.FLAV.IHEM) 


1 


BSHEM( 


:iBH. 


.FLAV45 , IHEM) 


1 


BSHEM( 


;iBH. 


.FLBOTH,IHEM) 


1 


BSHEM( 


:iBH. 


.MY, IHEM) 


1 


BSHEM( 


:iBH. 


.SHEM,IHEM) 


1 


BSHEM( 


;iBH. 


.SHEM2,IHEM) 


1 


BSHEM( 


:iBH. 


.QUAL.IHEM) 


1 


BSHEM( 


:iBH. 


.qUAL2,IHEM) 


1 


BSHEM( 


:iBH. 


.MVER.IHEM) 


1 


BSHEM( 


:iBH. 


.BPLUS.IHEM) 


1 


BSHEM( 


:iBH. 


.BZERO , IHEM) 


1 


BSHEM( 


;iBH. 


.PRPF,IHEM) 


1 


BSHEM( 


;iBH. 


.PRPB , IHEM) 


1 


BSHEM( 


:iBH. 


.PRSF.IHEM) 


1 


BSHEM( 


;iBH. 


.PRSB.IHEM) 


1 


BSHEM( 


:iBH. 


.DIPOL,IHEM) 


1 


BSHEM( 


:iBH. 


.HKST.IHEM) 


1 


BSHEM( 


;iBH. 


.NPALL,IHEM) 


1 


BSHEM( 


:iBH. 


.NPGOOD.IHEM) 


1 


BSHEM( 


:iBH. 


.NPCUT,IHEM) 


1 


BSHEM( 


;iBH. 


.NPGC , IHEM) 


1 


BSHEM( 


:iBH. 


.NKS , IHEM) 


1 


BSHEM( 


:iBH. 


.NWKSD.IHEM) 


1 


BSHEM( 


:iBH. 


.NWKSF.IHEM) 


1 


BSHEM( 


:iBH. 


.NLAMO , IHEM) 


1 


BSHEM( 


;iBH. 


.NWLAMBD,IHEM) 


1 


BSHEM( 


;iBH. 


.NWLAMBF.IHEM) 


1 


BSHEM( 


:iBH. 


.NNEUT,IHEM) 


1 


BSHEM( 


:iBH. 


.NPION,IHEM) 


1 


BSHEM( 


:iBH. 


.PICHRG.IHEM) 


1 


BSHEM( 


:iBH. 


.NKDECAY,IHEM) 


1 


BSHEM( 


;iBH. 


.NKFRAG,IHEM) 


1 


BSHEM( 


:iBH. 


.NWKD.IHEM) 


1 


BSHEM( 


:iBH. 


.NWKF.IHEM) 


1 


BSHEM( 


:iBH. 


.NPDECAY,IHEM) 


1 


BSHEM( 


;iBH. 


.NPFRAG.IHEM) 


1 


BSHEM( 


;iBH. 


.NWPROTD.IHEM) 


1 


BSHEM( 


:iBH. 


.NWPROTF,IHEM) 


1 


BSHEM( 


;iBH. 


.QKAMAX.IHEM) 


1 


BSHEM( 


:iBH. 


.PKAON , IHEM) 


1 


BSHEM( 


;iBH. 


.QKAON,IHEM) 


1 


BSHEM( 


:iBH. 


.PLEP , IHEM) 


1 


BSHEM( 


:iBH. 


.QLEP.IHEM) 


1 



Original b flavour NN 

improved b flavour NN (<0.5->b, >0.5->bbar) 
likelihood combination with opp. side info, 
rapidity mass 

scaled measured hemisphere energy- 
scaled measured hemisphere energy (2-BODY) 

hemisphere quality flag: 

N(rejected trk cuts)+100*N(interactions)+ 
1 , OOON (ambiguous hits) +100 , 000*N(ISRT=0) 
Continous hemisphere quality 0-(perf ect) lO(bad) 
Mass on secondary vertex 
B+ tagging 
BO tagging 

prim, vertex fragmentation Borisov prob 

prim, vertex B Borisov prob 

sec. vertex fragmentation Borisov prob 

sec. vertex B Borisov prob 

sec . vertex dipole moment 

Q*PMIN of tracks in B CMS 

# of particles 

# of particles past quality cuts 

# particles past IBP.NET > PNET.CUT 

# past IBP_NET>PNET_CUT + qual. requirement 

# recon. K_s with y>1.6 

sum of y-prob. funct for decay K_s 
sum of y-prob. funct for frag. K_s 

# recon. Lambda(uds) with y>1.6 

sum of y-prob. funct for decay lambdas 
sum of y-prob. functions for frag. lambdas 

# recon. neutrons with y>1.6 

# charged pions 

# pions charged weighted about SV 
Sum of kaon weights for trknet> 0.5 
Sum of kaon weights for < 0.5 

sum of y-prob. funct for decay K+/- 
sum of y-prob. functions for frag. K+/- 
Sum of proton weights for trknet > 0.5 
Sum of proton weights for trknet < 0.5 
sum of y-prob. functions for decay protons 
sum of y-prob. functions for frag, protons 
MAX. q*trknet in hemisphere for tight K. 
charge correl. coeff . best kaon cand. in hem 
charge of best kaon cand. in hemisphere 
charge correlation coeff . of best lep 
charge of best lepton cand. in hemisphere 
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* BSHEM(IBH_DEC3D,IHEM) 1 3D decay length 

* BSHEM(IBH_DEC3DS,IHEM) 1 error on 3 D decay length 

* BSHEM(IBH_BLIFE,IHEM) 1 B lifetime 

* BSHEM(IBH_BLIFES,IHEM) 1 error on B lifetime 

* BSHEM(IBH_CHRDIP,IHEM) 1 B-D charge dipole moment from GBDIPOLE 

* BSHEM(IBH_SVVAR1,IHEM) 1 Mean deviation of trks about SV point 

* BSHEM(IBH_SVVAR2,IHEM) 1 Mean square deviation 

* BSHEM(IBH_DIPPY,IHEM) 1 A dipole moment from the SVVAR-calc 
* 

* BSHEM(IBH_QRANK1,IHEM) 1 Charge of leading frag, track 

* BSHEM(IBH_QRANK2,IHEM) 1 Charge of 2ND rank frag track 

* BSHEM(IBH_qRANK3,IHEM) 1 Charge of 3RD rank frag track 

* BSHEM(IBH_YRANK1,IHEM) 1 Rapidity of leading frag track 

* BSHEM(IBH_YRANK2,IHEM) 1 Rapidity of 2ND rank frag track 

* BSHEM(IBH_YRANK3,IHEM) 1 Rapidity of 3RD rank frag track 

* BSHEM(IBH_KRANK1,IHEM) 1 Kaon net for leading rank frag track 

* BSHEM(IBH_KRANK2,IHEM) 1 Kaon net for 2ND rank frag track 

* BSHEM(IBH_KRANK3,IHEM) 1 Kaon net for 3RD rank frag track 

* BSHEM(IBH_PRANK1,IHEM) 1 Proton net for leading rank frag track 

* BSHEM(IBH_PRANK2,IHEM) 1 Proton net for 2ND rank frag track 

* BSHEM(IBH_PRANK3,IHEM) 1 Proton net for 3RD rank frag track 

* BSHEM(IBH_DELTA,IHEM) 1 delta - y-gap above/below trknet=0.5 

* BSHEM(IBH_QFSUM,IHEM) 1 Sum of frag, track charge 

* BSHEM(IBH_DFLBS,IHEM) 1 decay flavour net for B_s 

* BSHEM(IBH_DFLBO,IHEM) 1 decay flavour net for BO 

* BSHEM(IBH_DFLBP,IHEM) 1 decay flavour net for B+ 

* BSHEM(IBH_DFLLB,IHEM) 1 decay flavour net for LamB 

* BSHEM(IBH_FFLBS,IHEM) 1 fragmentation flavour net for B_s 

* BSHEM(IBH_FFLBO,IHEM) 1 fragmentation flavour net for BO 

* BSHEM(IBH_FFLBP,IHEM) 1 fragmentation flavour net for B+ 

* BSHEM(IBH_FFLLB,IHEM) 1 fragmentation flavour net for LamB 

* BSHEM(IBH_PRBBS,IHEM) 1 probability for B_s from all-flavour net 

* BSHEM(IBH_PRBBO,IHEM) 1 probability for BO from all-flavour net 

* BSHEM(IBH_PRBBP,IHEM) 1 probability for B+ from all-flavour net 

* BSHEM(IBH_PRBLB,IHEM) 1 probability for LamB from all-flavour net 

* BSHEM(IBH_PRBAIBS,IHEM) 1 B_s prob. from all-flavour net inc. opp. hem. info. 

* BSHEM(IBH_PRBAIBO,IHEM) 1 BO prob. from all-flavour net inc. opp. hem. info. 

* BSHEM(IBH_PRBAIBP,IHEM) 1 B+ prob. from all-flavour net inc. opp. hem. info. 

* BSHEM(IBH_PRBAILB,IHEM) 1 LamB prob. from all-flavour net inc. opp. hem. info. 

* BSHEM(IBH_DFLAV45,IHEM) 1 Optimal decay flavour net 

* BSHEM(IBH_LIKBS,IHEM) 1 Likelihood ratio prod. flav. tag for B_s 

* BSHEM(IBH_LIKBO,IHEM) 1 Likelihood ratio prod. flav. tag for B_d 
* 

* BSHEM(IBH_WEIGHT,IHEM) 1 SIMWAIT weight for physics/quality differences 

* between simulation and data 
* 

* BSHEM(IBHMC_BDEC2D,IHEM) 1 MC generated 2-d B decay length 
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BSHEMCIBHMC. 
BSHEMCIBHMC, 
BSHEMCIBHMC, 
BSHEMCIBHMC. 
BSHEMCIBHMC. 
BSHEMCIBHMC. 
BSHEMCIBHMC, 
BSHEMCIBHMC, 
BSHEMCIBHMC. 
BSHEMCIBHMC. 
BSHEMCIBHMC. 
BSHEMCIBHMC. 
BSHEMCIBHMC, 
BSHEMCIBHMC, 
BSHEMCIBHMC. 
BSHEMCIBHMC. 
BSHEMCIBHMC. 
BSHEMCIBHMC, 
BSHEMCIBHMC, 
BSHEMCIBHMC. 
BSHEMCIBHMC. 



.BDEC.IHEM) 1 

.DDEC , IHEM) 1 

.DPDECIHEM) 1 

.BV,IHEM) 3 

.DV,IHEM) 3 

.DPV.IHEM) 3 

.DST,IHEM) 1 

.BP, IHEM) 4 

.BKF , IHEM) 1 

.ND,IHEM) 1 

.OSCB.IHEM) 1 

.BLIFE.IHEM) 1 

.BMC, IHEM) 1 1 

.BMTOT,IHEM) 1 

.QTYP , IHEM) 1 

.PBKF.IHEM) 1 

.DRSKF,IHEM) 1 

.qCDRS , IHEM) 1 

.QNDRS , IHEM) 1 

.DRSP,IH) 1 

.SLD,IH) 1 



MC generated 3-d B decay length 
MC generated 3-d Dbar decay length 
MC generated 3-d D decay length 
MC generated B vertex 
MC generated Dbar vertex 
MC generated D vertex 
MC D* in hemisphere 
MC B 4-VECTOR 

MC kf code for weakly decaying B 
MC: # of weakly decaying D 
MC flag for oscillating B* 
MC B lifetime 

MC B charged decay multiplicity 

MC B charged+neutral decay multiplicity 

MC produced quark flavour 

MC kf code for primary B 

MC kf code for right sign D 

MC charged decay mult . of right sign D 

MC charged+neutral decay mult . of RSD 

MC P of the right sign D in the B CMS 

MC l=weak B has decayed to e or mu, O=not 



The track output array where IPART is the track number: 

Array Element Number of Description 

values 



BSPAR C IBP.HEM , IPART) 
BSPARCIBP.IYB, IPART) 
BSPAR C IBP.QUAL , IPART) 



BSPAR C IBP. 
BSPAR C IBP. 
BSPAR C IBP. 
BSPAR C IBP. 
BSPAR C IBP. 
BSPAR C IBP. 
BSPAR C IBP. 
BSPAR C IBP. 
BSPAR C IBP. 
BSPAR C IBP. 



.Y, IPART) 
PRP, IPART) 
PRS , IPART) 
BRP , IPART) 
.BZ, IPART) 
.B3D, IPART) 
BSRP, IPART) 
BSZ, IPART) 
,BS3D, IPART) 
.NET, IPART) 



BSPAR C IBP.DNET , IPART) 
BSPAR CIBP.NRP, IPART) 
BSPARCIBP.NZ, IPART) 
BSPAR C IBP.KST , IPART) 
BSPAR CIBP.THETST, IPART) 1 



hemisphere Cor jet) number 

index of the matching VECSUB track 
track quality: 20Cif isrt=0) + lOCif from 
interaction) + iCif has ambiguous hits) 
rapidity with respect to jet 
probability to fit primary vertex 
probability to fit secondary vertex 
Borisov probability in rphi 
Borisov probability in z 
Borisov probability in 3D 

Borisov probability in rphi w/r secondary vtx 
Borisov probability in z w/r secodary vtx 
Borisov probability in 3D w/r secodary vtx 
net output C — >0 for good fragmentation, 

— >1 for good B decay product) 
cascade D track net output 
number of rphi VD hits 
number of z VD hits 
momentum of track in B CMS 
helicity angle of the track 
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BSPARCIBP, 
BSPARCIBP 
BSPARCIBP 
BSPARCIBP, 
BSPARCIBP, 
BSPARCIBP, 
BSPARCIBP 
BSPARCIBP 
BSPARCIBP, 
BSPARCIBP, 
BSPARCIBP, 
BSPARCIBP, 
BSPARCIBP 
BSPARCIBP, 
BSPARCIBP, 
BSPARCIBP, 
BSPARCIBP, 
BSPARCIBP 
BSPARCIBP, 



.BSAU.IPART) 1 

.SV,IPART) 1 

.TRKE,IPART) 1 

.TRKM,IPART) 1 

.TRKP , IPART) 1 

.TRKQ , IPART) 1 

.TRKL, IPART) 1 

.ERRE, IPART) 1 

.DELPP, IPART) 1 

.DFLBS, IPART) 1 

.DFLBO, IPART) 1 

.DFLBP, IPART) 1 

.DFLLB , IPART) 1 

.FFLBS , IPART) 1 

.FFLBO , IPART) 1 

.FFLBP, IPART) 1 

.FFLLB, IPART) 1 

.MACK, IPART) 1 

.MACP , IPART) 1 



pointer to BSAURUS internal particle array 

1 — > trk in sec. vertex — > trk not in SV 

Track energy 

Track mass 

Track momentum 

Track charge 

Track length 

Error on track energy 

DeltaCP)/P 

charge/decay flavour correlation net for B_s 
charge/decay flavour correlation net for BO 
charge/decay flavour correlation net for B+ 
charge/decay flavour correlation net for LamB 
charge/f ragm. flavour correlation net for B_s 
charge/f ragm. flavour correlation net for BO 
charge/f ragm. flavour correlation net for B+ 
charge/f ragm. flavour correlation net for LamB 
MACRIB kaon net output 
MACRIB proton net output 



BSPAR C IBPMC.KF , IPART) 1 
BSPAR C IBPMC.KFP , IPART) 1 
BSPAR C IBPMC.TYP , IPART) 1 



MC Lund KF code 

MC Lund KF parent code 

MC code 

: unidentified 

1 : ordinary fragmentation 

2 : fragmentation partner of B-hadron 

3 : decay of primary Cnot weak) C 

4 : decay of D* 

5 : decay of weak B 

6 : decay of excited B 

7 : decay of weak C 

Ix : same as for x but decay of 

long living particles Ce.g. K_l) 
2x : same as for x but from 

hadronic interaction or photo conv. 
-/+ : b/bbar quark 

6 : content of primary B 

3,4,5,7 : content of weak B 



BSPAR C IBPMC.KST , IPART) 
BSPAR C IBPMC.THETST , IPART) 
BSPAR C IBPMC.TRKP , IPART) 
BSPAR C IBPMC.TRKE , IPART) 



1 MC momentum of track in MC B CMS 

1 MC helicity angle of the track 

1 truth track momentum 

1 truth track energy 



Other: 

Array Element 



Number of 



Description 
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values 



NTOT= total number of particles 
NCH= number of charged tracks 

Hemisphere status word ISAURUS 



ISAURUS 
ISAURUS 
ISAURUS 
ISAURUS 
ISAURUS 
ISAURUS 



(IHEM) 
(IHEM) 
(IHEM) 
(IHEM) 
(IHEM) 
(IHEM) 




1 
2 
3 
4 
5 



event not processed 

event accepted by user routine FORUSE 
rapidity algorithm successful 
rapidity energy > EYMIN (10 GeV) 
sec. vertex fit successful 
sec. vertex fit with convergence 
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